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The problem of membrane topology in the matrix model of M-theory is considered. 
The matrix regularization procedure, which makes a correspondence between finite- 
sized matrices and functions defined on a two-dimensional base space, is reexamined. 
It is found that the information of topology of the base space manifests itself in the 
eigenvalue distribution of a single matrix. The precise manner of the manifestation is 
described. The set of all eigenvalues can be decomposed into subsets whose members 
increase smoothly, provided that the fundamental approximations in matrix regular- 
ization hold well. Those subsets are termed as eigenvalue sequences. The eigenvalue 
sequences exhibit a branching phenomenon which reflects Morse-theoretic information 
of topology. 

Furthermore, exploiting the notion of eigenvalue sequences, a new correspondence 
rule between matrices and functions is constructed. The new rule identifies the matrix 
elements directly with Fourier components of the corresponding function, evaluated 
along certain orbits. The rule has semi-locality in the base space, so that it can be 
used for all membrane topologies in a unified way. A few numerical examples are 
studied, and consistency with previously known correspondence rules is discussed. 
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1 Introduction 



Topological properties of a system are often important in investigating the dynamics of 
the system. It seems certain that M-theory[l] has membranes as dynamical degrees of 
freedom. Furthermore, the only existing proposal for formulation of M-theory, namely the 
matrix model of M-theory [2], can be considered as an attempt to define quantum membrane 
theory [3, 4, 5]. More explicitly stated, it is a regularized version of membrane theory in 
lightcone gauge, dynamical variables becoming N x N matrices instead of functions defined 
on two-dimensional worldspace. 

But, at present, the topological properties of membranes in M-theory are not known. 
The concern of this paper is membrane topology 1 in the matrix model. It is believed that 
the matrix model can describe membranes of arbitrary topologies. However, there has been 
a problem: we do not know whether and how the information of the topology manifests 
itself in the matrix model. The cause for this problem lies in the manner in which the 
correspondence between matrices and functions has been given. There has been no unique 
rule that can deal with all membrane topologies. Instead, we have many different rules for 
different topologies [3, 4, 5, 6, 7], interrelationships between those rules being unclear. 

In this paper, we address this problem by reexamining the regularization procedure, the 
so-called matrix regularization. We shall show that the information indeed manifests itself in 
the eigenvalue distribution of a single matrix. The precise manner of the manifestation will 
be described. Moreover, we have constructed a new correspondence rule between functions 
and matrices which can be applied to all membrane topologies in a unified way. 

We start by discussing relevant aspects of matrix regularization, in section 2. We take 
the simple view that, classically, the matrix regularization is an approximation of continuum 
theory by a discretized theory. 2 This approximation between two theories is based solely 
on some fundamental large- N approximation formulae (2.7)-(2.9). They play a vital role 
in this paper. We also recall the well-known mathematical analogy between the matrix 
regularization and canonical quantization of systems with one degree of freedom, which will 
be our main tool in subsequent discussions. 

Then, in section 3, we turn to the investigation of membrane topology in the matrix 
model. Our basic observation is that, in order to study membrane topology, it suffices to 
consider the two-dimensional base space, which we shall term as the cr-space, not the shape 

x We use the word membrane topology to express the topology of a configuration of membranes in a time- 
slice. The topology is not of a single membrane but of a totality of membranes. Thus, membrane topology 
is classified by the numbers of membranes rn which has genus i = 0, 1, • • •. 

2 In quantum theory, at the same time, the matrix regularization is considered as a definition of continuum 
theory by a non-trivial limit of discrete theories. This is the reason why we should first treat the finite- N 
theory carefully. 
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of the membranes in the target space. This observation greatly simplifies the analysis, since 
it enables us to deal with only a single matrix, not many matrices. 

We base the discussion on the analog of the Bohr-Sommerfeld quantization condition. We 
shall show that, in the case where the fundamental approximations hold well, the eigenvalue 
distribution of a matrix has a particular structure. Namely, the set of all eigenvalues can be 
decomposed into subsets characterized by the following property: the eigenvalues in one of 
the subsets, when sorted, increase smoothly. We call these subsets as eigenvalue sequences. 
The grouping of the eigenvalues into sequences reveals a branching phenomenon of sequences. 
We find that the branching phenomenon, in turn, reflects certain Morse-theoretic information 
of topology of the cr-space. This is our answer to the above problem. Thus, the information of 
topology manifests itself, in the world of matrices, as a branching phenomenon of eigenvalue 
sequences. 

Furthermore, the notion of eigenvalue sequences enables us to construct a new correspon- 
dence rule between matrices and functions, which is the subject in section 4. The matrix 
elements are approximately equal to Fourier components of the corresponding function, cal- 
culated along appropriate orbits on the cr-space. The rule is analogous to the correspondence 
noticed by Heisenberg when he created Matrix Mechanics pursuing Bohr's correspondence 
principle [8]. There, the matrix elements of an observable, in the basis which makes the 
Hamiltonian diagonal, are equal to the classical Fourier components of the observable along 
the appropriate classical orbits on the phase space. We shall show that the fundamental 
approximation formulae hold well if the new correspondence rule holds. The correspondence 
rule contains the above-mentioned analog of the Bohr-Sommerfeld condition. This justifies 
the use of it in section 3. 

The new rule is semi-local in the cr-space, and consequently can be applied for all mem- 
brane topologies uniformly, in marked contrast with the previously known rules. This, in 
particular, enables one to construct functions corresponding to given matrices when the 
approximations are good. Using previous rules, one could only do the reverse, namely, to 
construct matrices corresponding to given functions. This is because one could not know 
the topology corresponding to the given matrices, and therefore could not choose the rule to 
be used. 

Apart from the unified treatment for all topologies, the new rule has the virtue that 
the identification of the matrix elements with Fourier components is direct, and so that the 
geometrical meanings of the matrix elements are clear. Our arguments are also relevant to 
the matrix model of type IIB string theory [9], since the same regularization is involved. 
Further, the same kind of mathematics as that of matrix regularization appears in such 
subjects as bound states of D-branes or non-commutative field theory. Ideas in this paper 
may find some applications in those subjects. 

A few illustrative numerical examples are given in section 5. The consistency between 
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our new rule and previous rules is checked by studying them. Finally, we conclude with some 
discussions in section 6. 



2 Matrix regularization 

Let us briefly recall the matrix regularization procedure from our viewpoint. Although it 
is supermembrane theory in eleven dimension [10] that is relevant to M-theory, we consider 
bosonic membrane theory for simplicity of presentation. 

Firstly, we shall describe the continuum theory. We parametrize the membranes by three 
parameters (r, a 1 , a 2 ) = (r, cr). Then, the geometrical shape of membranes in spacetime is 
described by the coordinate functions rr M (r, cr). In lightcone gauge formalism, r is chosen to 
be equal to x + , and cr is chosen so as to make the area of a domain in the cr-space proportional 
to total p + contained in the domain. Here, we denote the momentum density vector of the 
membranes by p^. The canonical variables of the system are transverse coordinates and 
momenta (which are functions defined on the cr-space) as well as zero modes, 

x a (<r),^(<r);X-,-P + . (2.1) 

The Hamiltonian is given by 

S=M /W« A , (2 . 2) 

with Lie brackets 

df dg df dg 

if, 9} 



da 1 da 2 da 2 da 1 ' 

where / and g are functions on the cr-space. We have also introduced a conventional constant 
[cr] which is the total area of the cr-space, [cr] = / d 2 cr. The remaining ingredients of the 
theory are the phase space constraints 

{x a ,p a }(a) = 0, (2.3) 

and its global version. They correspond to the local symmetry of the lightcone gauge theory 
under reparametrization by area-preserving diffeomorphism (APD) on the cr-space. This is 
a local symmetry, because one can perform reparametrization by different APD for different 

T. 

Secondly, we shall give the regularized theory. The canonical variables are NxN matrices 
as well as zero modes, 

x a ,p a ;X~,-P + . (2.4) 



3 



The Hamiltonian is given by 

H = NTi— ^ — — — — , 2.5 

and the constraints are, 

[x a ,p a ] = 0, (2.6) 

where [ , ] is a commutator of matrices. 

Now, we turn to the explanation of the matrix regularization. The following fact is essen- 
tial: there exists a correspondence between appropriate functions on the cr-space f(a),g(cr), ■ ■ 
and matrices /,#,••• such that the fundamental approximation formulae 

±ff(*)** « iTr/ (2.7) 

fg « fg (2.8) 
{Tg} « -^[/^] (2.9) 

hold. Here, we denote by {/, (?} and the matrices which correspond to the functions 
{/, g}(<r) and f(cr)g(cr) ) respectively. 3 The larger is N, the better is the approximation. 
From these formulae it follows that the continuum theory, defined by (2.1)-(2.3) can be 
approximated by a regularized theory defined by (2.4)-(2.6). We stress the importance of 
above formulae. They are almost the definition of the matrix regularization. 

Since Lie brackets and matrix commutators both obey the Jacobi identity and antisym- 
metry, the important advantage of matrix regularization follows. Namely, the regularized 
theory has local symmetry under the transformation 

x a '( T ) = [/(^^(^[/(r)- 1 , p a '{r) = U{T)p a (T)U(T)-\ (2.10) 

where U{r) is an arbitrary matrix which is a function of r, corresponding to the APD 
symmetry in continuum theory. 

The matrix regularization procedure is analogous to the quantization of a system which 
has one degree of freedom, as is well known. The analogy can be summarized as, 



Canonical quantization 


Matrix regularization 


(x,p) 


(<rV 2 ) 


Canonical transformation 


Area-preserving diffeomorphism 


{ , }P.B. - , ] 




h 


m 

2nN 



(2.11) 



3 Maybe we should add the linearity of the correspondence, / + g = f + g, for the sake of completeness. 
We have omitted it since it holds trivially in all our discussions. 
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where { , }p.b. is the usual Poisson brackets. We shall motivate our discussion by this 
analogy in section 3. 

We conclude this section with discussions on the previously known correspondence rules. 
We first recall the general manner the rules are formulated. We must, first of all, fix topology 
of the cr-space. After that, we consider a basis in the vector space of all functions defined 
on the cr-space. Then, we define an appropriate basis in the vector space of all N x N ma- 
trices, and postulate a correspondence between it and the basis in the space of the functions 
appropriately truncated. The rules are, finally, justified by checking that the fundamental 
approximations (2.7)-(2.9) hold well for large N by them. 

This manner has made difficult to consider whether and how membrane topology mani- 
fests itself in the matrix model. In particular, one can expand an arbitrary matrix by basis 
referring to any particular topology. This fact, at first sight, seems to suggest that a config- 
uration of matrix model could be interpreted as membranes of arbitrary topology, and there 
would be, therefore, no information of topology in the matrix model. 

This is not necessarily true. Even if one can formally expand some matrices by a basis 
referring to a particular topology, the fundamental approximations may not work at all. 4 In 
our perspective, that the matrix regularization is an approximation scheme, we cannot, then, 
interpret the matrices as membranes of the particular topology. Information of topology 
may be hidden in the matrices in this way. Through sections 3 and 4, we shall see indeed 
that, provided that the approximations are good, the information reflects in the eigenvalue 
distribution. 

3 Membrane topology and matrix regularization 

In this section, we show that the information of membrane topology manifests itself in the 
matrix model. Before explicit description of the manner of the manifestation, let us give 
some basic observations. 

If one wishes to specify the complete shape of membranes in the target space, one needs 

4 This may be expected from the previously known rules. Let us, for example, imagine a smooth function 
defined on a torus. One can construct the corresponding matrix using the basis for torus topology. One can 
then expand the matrix by the basis (in the space of matrices) corresponding to the topology of a sphere, and 
construct a function defined on a sphere. We expect that the resulting function would have discontinuity or, 
in any case, some singularity (see subsection II. C of [11]). This implies that the function varies considerably 
in a very small length scale. Therefore the approximations may well be no good, since, in general, the smaller 
the length scale of the variation of functions, the larger must be N in order that the approximations are 
good. 

However, it is difficult to characterize precisely, using only previously known correspondence rules, when 
the approximations break down. Hence, it has not been clear if this picture is indeed right. 
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information of many (that is, roughly speaking, as many as the dimension of the target 
space) functions. This would imply that one should study many matrices in the matrix 
model. However, the information of membrane topology, or at least the information of 
topology of the cr-space, is contained in one generic function defined on the cr-space, as 
is strongly suggested by Morse theory. We choose, as our basic strategy, to consider the 
latter information. Then, we shall seek in a single matrix the information of topology of the 
cr-space. 

There is another point we would like to discuss. It is most natural to identify functions 
which are transformed into each other by APD transformations. We shall identify those 
matrices which are transformed into each other by similarity transformations, since (2.9) 
tells us that the counterpart of the APD transformation is the similarity transformation. 
This has some non-triviality, since it may happen that the identification is only allowed 
approximately. Nonetheless, we shall carry out the identification, because that the APD 
symmetry survives as (2.10) is the most important advantage of the matrix regularization. 
This identification and our strategy, to consider the topology of the cr-space, act together 
to greatly simplify the analysis. Since one can always diagonalize a single matrix, we can 
concentrate on the eigenvalue distribution of the matrix. 

Having explained our basic strategy, we shall now proceed to investigate the manner of 
the manifestation of membrane topology in the eigenvalue distribution. 

First, let us consider how one can read off the information of topology from a function in 
an APD invariant way. We choose an arbitrary generic function f(cr). 5 It could be one of 
the transverse coordinates, for instance. The function is fixed, throughout our discussion, as 
a kind of reference. Thus, we shall use /, shortly below, as both an analog of the Hamiltonian 
in canonical quantization and a Morse function. As a natural APD-invariant concept with 
a given function /(<x), we introduce an ordinary differential equation (ODE) 

jf = <"•»■ < 31 > 

drawing analogy to the Hamiltonian equation of motion with a given Hamiltonian function 
H(x, p), which is invariant under canonical transformations. This ODE governs the motion of 
points of the cr-space. 6 Thus, we envisage an auxiliary Hamiltonian-like dynamical system 
with the cr-space as its phase space and with / as its Hamiltonian. Since / is conserved 
along the motion by the identity {/, /} = 0, an orbit of this equation is a part of an equal-/ 
line in the cr-space. It will form a closed loop because of the compactness of the cr-space. 

5 We use the word generic in the sense of Morse theory: we avoid degenerate functions, constant functions 
for instance, which can be changed into generic functions by arbitrarily small perturbations. 

6 We note that the independent variable of the ODE, t, is just a mathematical tool to substantiate the 
analogy to canonical quantization. It has nothing to do with physical time coordinates of membrane theory. 
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Figure 1: The cr-space which has topology of a torus. The height in the figure is the reference 
function. Some orbits of the ODE der/dt = {cr, /} are drawn, which form closed loops that 
are, in turn, parts of the equal-/ lines. If one gradually increases the value of /, one observes 
a branching phenomenon of the orbits: appearing, branching, merging, and disappearing 
processes at the points A, B, C, and D, respectively. Depicted is essentially the cr-space, so 
that the horizontal directions of the figure have rather arbitrary meanings. If one wishes, 
one can also give definite meanings to the horizontal directions by interpreting this figure as 
the geometrical shape of a membrane in the target space, and the reference function as one 
of the coordinate functions. 



Then, if we scan the cr-space by gradually increasing the value of /, we will observe 
branching processes of these orbits. There are four types of these branching processes: 
appearing, disappearing, branching and merging. Let us consider, for a typical example, 
the situation depicted in Fig. 1. The membrane topology is that of a torus. The reference 
function / is chosen to be the height in the figure, and some orbits of (3.1) are drawn. In this 
example, at the points A, B, C, D, the orbits appear, branch, merge, disappear, respectively. 
We can read off the information of topology from these processes. This is just the well- 
known idea of Morse theory. In particular, we obtain the Euler number of the cr-space, by 
subtracting the total number of the branching and merging processes from the total number 
of the appearing and disappearing processes. 

Now, we shall show that this analysis of topology in the world of functions has a coun- 
terpart in the world of matrices. The analogy of the matrix regularization to canonical 
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quantization is useful here. In the latter, the Bohr-Sommerfeld quantization condition de- 
termines the eigenvalues of the Hamiltonian operator H from the classical Hamiltonian 
function H(x,p) defined on the phase space (x,p). Namely, we draw classical orbits in (x,p) 
space, that are parts of equal- H lines, so that the areas of the domains between two adjacent 
orbits are equal to 2irk. Then eigenvalues of H are given by the values of H at these orbits. 
Here, we shall exploit the analogy, which is summarized in (2.11), and state the analog of 
the Bohr-Sommerfeld condition. Namely, we draw orbits of (3.1) in the cr-space so that the 
areas of the domains 7 between two adjacent orbits are equal to [cr]/N. Since [a] is the total 
area of the cr-space, this simply means that we divide the cr-space into N parts of equal 
area. Eigenvalues of / are then given by the values of / at these orbits. We assume this rule 
to hold. We shall justify the assumption in section 4. 

If we apply this rule to the case in Fig. 1, then the eigenvalues of / can be grouped into 
four subsets each of which corresponds to the family of the orbits belonging to (i) the region 
from the point A to the point B, (ii) the left branch of the torus from the point B to the 
point C, (iii) the right branch of the torus from the point B to the point C, (iv) the region 
from the point C to the point D, respectively. We call these subsets as eigenvalue sequences. 
For large enough N, eigenvalues belonging to each sequence have the following property. If 
we sort the eigenvalues contained in a sequence in increasing order of their values, and make 
a graph plotting the values of them versus their order, then the plotted points can be linearly 
approximated locally. To put it short, the eigenvalues in a sequence increase smoothly. It 
should be clear that, in general, if we do not group the eigenvalues properly, then the graph 
become zigzag-shaped and the above property is lost. In section 4, we see that this linear 
approximation is essential in order the fundamental approximations (2.7)-(2.9) to hold. 

The eigenvalue sequences should exhibit the same branching phenomenon as that of the 
orbits. For the example of Fig. 1, the sequence (i) appears and then branches into the 
sequences (ii) and (iii). They merge into the sequence (iv), and finally (iv) disappears. 
It is clear that all these considerations work the same in general cases other than that 
of Fig. 1. Thus, the information of membrane topology manifests itself in the branching 
phenomenon (which consists of appearing, branching, merging, disappearing processes) of 
eigenvalue sequences. A few examples, including the case similar to the situation in Fig. 1, 
are given in section 5. 

7 The area of a domain means here area in the cr-space not in the target space. Its physical meaning is 
the total p + contained in the domain, apart from a conventional factor, by the gauge choice made in the 
lightcone gauge formalism. 
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4 The new correspondence rule 



In this section, we present a new correspondence rule between matrices and functions, and 
then show that the fundamental approximations (2.7)-(2.9) stem from the rule. 

We choose an arbitrary generic function / and fix it as a reference, as in section 3. The rule 
is formulated in such a way that the representation of matrices is so chosen that the matrix 
/, corresponding to the function /, is diagonal. For simplicity of notation, we shall consider 
the case where only one eigenvalue sequence is present. We explain the generalization later 
in this section. 

We first give the rule to determine the diagonal matrix /. To this end, we set up some 
notations. We again consider ODE (3.1) 

A solution of this ODE is periodic, the point of the cr-space circulating on a loop which is 
part of an equal-/ line. We shall denote its period, as a function of /, by T(f). We sort the 
eigenvalues of / in increasing order, and call them f n , 

■■■<fn-l< In </„+!<•••■ (4.1) 

To be specific, we choose the representation such that 

/ = diag(- • • , / n _i, /„, f n+1 , ■ ■ •)• (4.2) 

The relation between the function / and the matrix elements f n is the analog of the Bohr- 
Sommerfeld quantization condition stated in section 3. If N is sufficiently large, the rule can 
be formulated as, 

L-fn~(m- n)Q _ f l , (4.3) 

2 



TV 7^/m±Jk) 



when \m — n\ is small. We have used that for two nearby loops, one at / and the other at 
f + Sf, the area 5S between them can be approximated by 
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We can construct f m satisfying (4.3) directly by the following method. We first define 
S(f) = j (1/T(f))df. The value of S(f) runs from to [<r] in this case where there is only 
one eigenvalue sequence. We then consider the inverse function f(S), and set f m = f(S m ), 
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where S m are determined by S m+ i — S m = [<r]/N up to a constant shift. The shift should 
be of order 1/JV for consistency. 8 

Having stated the correspondence rule for the reference function, we next turn to the 
correspondence rule for an arbitrary function g. We denote the matrix elements of the 
corresponding matrix g by g mn . When \m — n\ is small, g mn is equal to the Fourier component 
of order m — n of the function g{cr{t)). Here, <r(t) denotes the solution of (3.1) along which 
the function /(<x) takes the (constant) value (f m + / n )/2. To obtain explicit formulae, we 
define the Fourier components g s (f) by 

g(tr(t)) = £ 9s (f) e W)«, (4.5) 

where the parameter / denotes the value of the function f(cr) along the solution cr(t). We 
then set, 9 

/ fm "I" fn\ f a r\ 

gmn = g m -n[ " • (4.6) 



2 

We also require that when \m— n\ gets larger, the value of g mn falls off rapidly. This condition 
naturally conforms with (4.6), provided that the function g is sufficiently smooth and N is 
sufficiently large. 

A comment to the rule (4.6) is in order. We have freedom to change the orbit cr(t) by 
translation of t. The amount of translation is a function of /, which we denote by At(f). 
By this transformation, g s (f) becomes 

Therefore, g mn changes into 

c*<™-">* 9mn . (4.8) 

This freedom has a counterpart in the world of matrices. Namely, we can change m-th 
eigenvector by a phase factor e lSm . By this transformation, g mn becomes 

e ^-^)g mn . (4.9) 



8 We can determine the shift by setting S\ = [<r]/(2N) for an eigenvalue sequence beginning with an 
appearing process. This is analogous to the 1/2 in the Bohr-Sommerfeld condition § pdq = (n + (l/2))2nh. 
The justification for the above rule comes from the fact that (4.11) holds at one more higher order in 1/N 
by this rule. In other words, the rule is just the midpoint rule for numerical integration. Similar rule exists 
for an eigenvalue sequence ending with a disappearing process. 

9 This relation is the direct analog of the correspondence, in semi-classical region, between quantum matrix 
elements and Fourier components along classical orbits, first introduced in [8]. Also, formulae which bear 
some resemblance to ours appear in [12], where a correspondence between membrane theory and Matrix 
String Theory is considered. See also [13]. 
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Comparing (4.8) and (4.9), we find that if 

S n - 5 rn « (m - n) ^ At(^±A) (4.10) 

holds, then the two transformations are approximately identical. We can construct 5 m satis- 
fying (4.10) from given At(f), provided that At(/) is sufficiently smooth and N is sufficiently 
large. 

Equations (4.3) and (4.6) constitute, then, our new correspondence rule. We shall now 
deduce the fundamental approximations (2.7)-(2.9) from the new rule. 

By (4.3), we have divided the cr-space into N domains around orbits along which f(cr) 
takes the values /i, • • • , /;v. We can evaluate the integral of an arbitrary function / g(er)d 2 er 
approximately, by summing up the average values of g(tr) on these loops multiplied by the 
areas of each domains. Since (4.5) and (4.6) tell us that the average value of g on the orbit 
along which f(cr(t)) = f n is g nn , and since each area is equal to [cr]/N, we obtain 

J g(a) d 2 <r*Y.9nr)jj-, (4.11) 

which is nothing but (2.7). 

We next consider multiplication of matrices constructed by (4.6) 

(gh) mn = 9ml hn = 9m-l (~ ^ ~) kl ~ n (~ \ ~) " ( 4 ' 12 ) 

Since g m -i and g\- n fall off rapidly when \m — l\ and |Z — n| are large, respectively, the terms 
in which / is not far away from m or n dominate the summation. Then, by (4.3), neglecting 
higher order terms in 1/N, we can replace both (f m + f t )/2 and (ft + f n )/2 by (f m + f n )/2. 
We have, therefore, 

(gh)mn ~ J2 9m-l (~ | ~ ) fef-n (~ | ~ ) = Wm-n ^^y^ ) = ^m„, (4-13) 

where the second equality is the convolution law of Fourier series. Thus, the matrix corre- 
sponding to the multiplication of the two functions approximately coincides with the multi- 
plication of matrices corresponding to the functions. We have derived (2.8). 

We have just seen that, to the leading order, the multiplication of the matrices is commu- 
tative, since gh = hg, as a matter of course. Incorporating one more higher order terms in 
1/N, we shall evaluate the non-commutativity of the matrices. Thus, from (4.12), we have 

(fh) mn - E [9rn-l{^^) + (/*«-» (^^) + ^=^-n) ■ (4-14) 
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Here, we set g f s (f) = dg s /df. 10 We have omitted the value of / at which g' or h! is evaluated, 
since that does not affect the results to the order we are working. By (4.3), it follows that 



N T\2 \ \ 2 NT, 
where we have introduced new dummy indices u — m — l,v — I — n. Then, finally, we have 

([g, h]) mn ~ ( E + {vg' u )hv) -{g^ hu 

JV 1 Z \u+v=m-n / 

= [ 4^(-^9u)K + g' u (vh v )). (4.15) 

In order to understand the relation of the last expression to the function {g, h}(cr), it is 
instructive to consider the special case h — f. Namely, we consider the case in which one of 
the functions is the reference function. In that case we have by (4.3), 

[§, f]mn = gmn(fn ~ fm) ~ ^^( n ~~ m )9m-n, (4.16) 

which is the special case of (4.15). On the other hands, the Lie brackets between g and / 
can be expressed by a solution of (3.1) as, 



(4.17) 

<X(t)=<7 



where the total derivative with respect to t is taken at the point where the Lie bracket is 
calculated. Then, from the definition of the Fourier component g s , (4.5), we get 

- — - 2n 

{9,f} mn = i ( m - n )Y 9m - n - ( 4 - 18 ) 



Comparing with (4.16), we obtain 



{9j}mn ~ -i^Tn-ig, fUn, (4.19) 



the special case of (2.9). 



10 We choose the orbits in (4.5) smoothly, so that dg s /df is well-defined. 
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The last expression in (4.15) and the above derivation of (4.19) suggest the natural 
generalization. We reinterpret the independent parameter of the ODE, t, as a function 
defined locally on the cr-space. Then, from (4.17), we find 

{«,/} = f = l, (4-20) 

which means that we can consider that (t, f) as canonically conjugate variables in terms of 
the analogous canonical formalism. It follows that, 

The definition of the Fourier components g s , (4.5), is now interpreted as the representation 
of g as a function of (t, /) 

+ 00 / 2tt \ 

g(tj)= J2 9sUy { ~ )S) \ (4-22) 

s=— oo 

Substituting (4.22) and the similar formula for h into (4.21) we get, 11 

({g, h}) s = £ (i^ug u )h' v - g' u (i^vh v ). (4.23) 

u+v=s 

(Terms in which (d/df) t acts on 1/T(f) cancel out.) By comparing this expression with 
(4.15), we finally prove (2.9), 

{9, h} mn « -i-^-id, h]mn- (4.24) 

Up to this point, our derivation has been confined to the case where there is only one 
eigenvalue sequence. The extension to the general case where there are several eigenvalue 
sequences is easy. Namely, we apply (4.3) and (4.6) within each sequences separately. They 
determine the matrix elements between eigenvectors belonging to the same sequence. We 
then set remaining matrix elements, that is, matrix elements between eigenvectors which 
belong to different sequences, to zero. Above derivations of the fundamental approximations 
work just the same. 

This argument means that we can concentrate on the behaviour of functions on one 
branch of the cr-space, ignoring the behaviour on other branches. Also, since in our argu- 
ments the matrix element g mn falls off rapidly when \m — n\ gets larger, we can ignore the 

11 Technically, that t is defined only locally poses a problem. However, we can cope with it easily by 
introducing patches on each of which t is well-defined, and considering the relation between the patches. 
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behaviour of the functions at the place differing much in the value of /. These properties 
render our new rule a semi-local nature. That is, both the rule and the approximations work 
locally in the direction / changes. This situation is somewhat reminiscent of the uncertainty 
principle in the analogous quantum mechanical case. We have chosen the representation to 
make / diagonal. This choice achieves minimum uncertainty in /, and, at the same time, 
makes the conjugate variable t maximally uncertain. That our rule can be applied to any 
topology may be considered as a direct consequence of this semi-locality. 

The linear approximation (4.3) has been essential in the machinery of the derivations of 
(2.7)-(2.9). Therefore, it seems that the linear approximation, hence the existence of the 
eigenvalue sequences is necessary in order that the approximations are good. Also the use 
of the analog of the Bohr-Sommerfeld condition in section 3 is justified, since the condition 
is nothing but (4.3). 

Unfortunately, it seems that our new rule does not apply to the following exceptional 
quantities: matrix elements near branching and merging processes. Our rule is essentially 
a WKB approximation. In the immediate vicinity of the branching and merging processes, 
there should be tunneling effects which make the WKB approximation unreliable. Consider 
an analog problem in quantum mechanics, that is, the motion of a particle in the double-well 
potential. It is possible to deal with each well separately semi-classically, for sufficiently small 
H, and for generic energy levels. Indeed, tunneling amplitudes between the wells in general 
are negligibly small, behaving like exp (— 0(l)/h). However, for those rare energy levels 
which have energy close to the value of the potential at the local maximum, the tunneling 
amplitudes are not negligible. The break down of our rule could also be expected from a more 
direct argument. The solution of (3.1), <r(t), in the vicinity of the branching and merging 
processes, spends most of the time near the branching point, moving very slowly. Then, even 
if g(cr) is a smooth function, g(<r(t)) might develop singularity. Then, the validity of the 
condition used in our argument, that g mn is negligible for large \m — n\, might be questioned. 

5 Examples 

In this section we shall present three examples. In the first example, by an analytical calcu- 
lation, we show the equivalence between our new rule and the previously known rules. Both 
diagonal and off-diagonal matrix elements are compared. In the remaining two examples, 
our purpose is mainly to illustrate the notion of eigenvalue sequences. We calculate numer- 
ically eigenvalues of matrices constructed by the previously known rules. We represent the 
resulting eigenvalue distribution in a method such that the structure discussed in section 3, 
namely the eigenvalue sequences and their branching phenomenon, can be easily seen. We 
confirm that the branching phenomenon of the eigenvalue sequences coincides with that of 
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Figure 2: The cr-space of spherical topol- 
ogy. The height is the reference function. 
The orbits of (3.1) appear at A and dis- 
appear at B. 
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Figure 3: Plot of the eigenvalues and the 
difference of the eigenvalues of the matrix 
corresponding to the height in Fig. 2. The 
eigenvalue sequence appears at A and dis- 
appears at B. 



the orbits of the ODE (3.1). We further numerically compute the eigenvalues by our rule 
(4.3), and compare them with those calculated by the previously known rules. 

Example 1 We consider the cr-space which has topology of a sphere. We represent the 
cr-space as an unit sphere in £, rj, C-space, 

e+v 2 +c = i (s.i) 

with the area element given by 

dS = sin 9dBd<j>, (5.2) 

where 9 and <fi are polar coordinates defined by ( = cos 9, £ = sin 6? cos 0, 77 = sin # sin 0. 
Then the Lie brackets are {£,77} = C> ' ' '■ We choose the simple reference function f — C- 
Fig. 2 represents the cr-space and the reference function. The orbits of (3.1) appear at the 
point A and disappear at the point B. 

We first construct the matrix ( corresponding to the function (, by our new rule. The 
area of the domain ( < (' is given by 

^4tt. (5.3) 
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Then, by the analog of the Bohr-Sommerfeld quantization condition, or (4.3), we obtain 12 

(Ci,---,a) = (-l + ^,-l + ^,-"-,l-^). (5.4) 

We further construct the matrices £ and fj, corresponding to the functions £ and rj. The 
solutions to the ODE (3.1) can be explicitly written as, 

+ = V 1 " C 2 ^ (5.5) 

(e-«7)(f) = y/T^?e-*. 

Then, from (4.6), the only non-zero matrix elements is, 



(i + if)) 



m+l,m 



2 



1 - (^^Ild^ j = y'l _ _JL_ - y ) = (£ - ^) m>m+1 . (5.6) 



We shall now compare these results with those obtained from the previously known rules. 
The rule for the spherical topology reads [3, 4, 5], 



^ywzi**' ^Va^tt^ ^Viv 2 ^ 2 ' (5 ' 7) 

where l x ,l y ,l z are generators of the representation of SU(2) with spin I = (N — l)/2. Since 
eigenvalues of l z are {—I, —/ + !,•••,/}, we have, 



(Ci,---,6v) = 



N-l 4 N + l 4 AT — 1 



N 2 -l 2 ' V N 2 - 1 2 ' ' V N 2 



(5.8) 



which coincides, for large N, with the result of our new rule, (5.4). Further, it is well known 
that in the basis where l z is diagonalized, l x and l y have matrix elements only between the 
eigenvectors corresponding to adjacent eigenvalues. The expression for the non-zero matrix 
elements are, 



< l' z + l\(l x + il y )\l' z >= y/l(l + 1) - I'M + 1) =< Q(l x - il y )\l' z + 1 > (5.9) 

where we have denoted by \l' z > the eigenvectors of l z belonging to the eigenvalue l' z . Thus, 
the result of the previously known rule is, 



+'»)), n + l.n. : I/' ~J^—~ {■"' W~) ( m ~— 7T~ ) - (L - "lirn.m + L- (5.10) 



12 See also footnote 8. 
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Figure 4: The cr-space of spherical topology with a more interesting branching phenomenon 
of orbits than Fig. 2. Appearing, branching, (first) disappearing and (final) disappearing 
processes occur at the points A, B, C and D, respectively. 

These matrix elements are also approximately equal to (5.6). The agreements of our rule with 
the previously known rule for the simple functions £, f], ( imply agreements for more general 
functions which can be constructed by multiplying £, 77, £ finite (much less than N) times. 
The reason for this is that the approximate equality between multiplication of functions and 
that of matrices, (2.8), is valid for both rules. 13 

The eigenvalues Q and the difference Cm — of the eigenvalues are given in Fig. 3. We 
see that the eigenvalues consist of one eigenvalue sequence. The sequence appears at the 
point A and disappears at the point B in Fig. 3. They correspond to the branching points 
of orbits A, B in Fig. 2. 

Example 2 We treat another case of spherical topology, which exhibits a more interesting 
branching phenomenon of eigenvalue sequences than the previous example. Perturbing the 
reference function considered there, we here consider the reference function of the form 



The reference function and the cr-space are schematically depicted in Fig. 4. Orbits of (3.1) 
appear at the point A and then branch into two families at the point B. Then, the orbits 

13 For torus topology, similar argument as in this example, using simple functions such as cos er 1 or sine 2 
(see example 3 for definitions), has a tricky aspect since these simple functions are degenerate functions in 
the sense of Morse theory. 



f(a) = a( + bt + c£ 2 . 



(5.11) 
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Figure 5: The plot of the eigenvalues and their difference, of the matrix corresponding 
to the reference function (5.11) given in Fig. 4. The eigenvalue distribution is calculated 
both by the previously known rule (open squares), and our new rule, namely, the analog 
of the Bohr-Sommerfeld condition (crosses). They agree almost completely. The branching 
phenomenon for eigenvalue sequences is the same as that of the orbits in Fig. 4. Horizontal 
lines signify critical values of / at which the processes in the branching phenomenon take 
place, calculated directly from /. 
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belonging to the right branch disappear at the point C, and finally the orbits belonging to 
the left branch disappear at the point D. 
The corresponding matrix / is given by 

if one uses the previously known correspondence rule (5.7). We have computed numerically 
its eigenvalues, in the case a = 1,6 = 2, c = 6, with N = 40. We have also obtained the 
eigenvalue distribution from our new rule. To this end, we have computed numerically the 
area of the cr-space as a function of the height / for each branches of the cr-space. Then, by 
the analog of the Bohr-Sommerfeld condition (4.3), 14 we have calculated the eigenvalues of 
/• 

We represent the eigenvalues by the following method to see the information of membrane 
topology. Firstly, we sort the eigenvalues in increasing order, 

Ai < A 2 < ••• < Ajv. (5.13) 

They are given in Fig. 5. In order to see the branching phenomenon clearly, it is useful to 
plot also the difference of the eigenvalues Aj+i — Aj. By the plot one finds that from the point 
B to the point C, the plot of Aj is zig-zag shaped. Thus, the plot of Aj gives a juxtaposition 
of four eigenvalue sequences. We see the same branching phenomenon of the sequences as 
that of the orbits in Fig. 4. The agreement between our new rule and the previously known 
rule is remarkable. 

Example 3 We consider the cr-space which has topology of a torus. The cr-space can be 
represented by [0, 2n) x [0,2n), where periodic boundary conditions are understood. We 
choose the reference function to be 

/(cr) = a cos a 1 + 6 cos a 2 . (5-14) 

We assume that a ^ 6, a ^ 0, 6 ^ 0, in order to avoid degenerate reference functions. To 
be specific we choose < a < 6. The reference function is represented in Fig. 6. It has 
essentially the same feature as the reference function in Fig. 1. At the points A, B, C, D the 
function / takes the critical values —a — 6, a — 6, —a + 6, a + 6, respectively. 

In the previously known correspondence rule for torus topology, one postulates [6], 

e^ 1 = h u e ia2 = h 2 , (5.15) 

14 See also footnote 8. 
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Figure 6: (a) Contour plot of the reference function / = a cos a 1 + b cose 2 , with < a < b. 
(b) Schematic picture of the cr-space and the reference function. The contours, i.e., the orbits 
of (3.1) appear, branch, merge and disappear at the points A, B, C, and D, respectively. In 
(a) the orbits are so written that the areas of the domains between two adjacent orbits are 
W/N. 

where hi and h 2 are the well-known N x N matrices which satisfy the relation h\h 2 = 
h 2 h 1 exp (i2n/N). Then, it follows that 

f = ^(h 1 + h\) + b -(h 2 + hl). (5.16) 

We have computed the eigenvalue distribution of this matrix numerically, in the case a = 
1,6 = 3, with N = 30. We can also calculate them by the new rule as we have done in the 
previous example. 15 The results by the two methods are given in Fig. 7. They agree well, 
except at the vicinity of the branching process at the point B. The reason for the discrepancy 
is noted at the end of the previous section: our new rule should not be trusted in the vicinity 
of branching processes. We can trust the previously known rule, on the other hand, since 
the fundamental approximations (2.7)-(2.9) are guaranteed by the rule (5.15), irrespectively 
of the branching phenomenon. 

In Fig. 7, we see a sequence, which appears at the point A and branches into two sequences 
at the point B. Then, the two sequences merge at the point C, and finally the last sequence 

15 Due to the poor knowledge in the vicinity of the merging and branching processes discussed at the end of 
section 4, we have two (or rather one due to the symmetry of the present example) undetermined parameters 
of order 1/N mentioned in footnote 8. We have fixed the order \/N parameter by comparison to the result 
of the previously known rule. 
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Figure 7: Same as Fig. 5 but for the reference function (5.14) given in Fig. 6. The topology of 
the cr-space is that of a torus. The eigenvalues calculated by the new rule and the previously 
known rule agree well except at the immediate vicinity of the point B or the point C. The 
branching phenomenon of the eigenvalue sequences is the same as that of the orbits in 
Fig. 6. The (approximate) degeneracy of eigenvalues from the point B to the point C is only 
accidental, being result of the symmetry of the /. 
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disappears at the point D. These branching processes of sequences directly correspond to 
the branching processes of the orbits A, B, C, D in Fig. 6. 

6 Conclusions and discussions 

In this paper, we have clarified some elementary but unknown features of the matrix regu- 
larization procedure. We have worked under the simple view that it is an approximation of 
a continuum theory by a discrete theory. The approximation between two theories is based 
solely on the fundamental approximation formulae (2.7)-(2.9). We have constructed a new 
geometrical correspondence rule between matrices and functions. We have shown the valid- 
ity of the rule directly by deriving the fundamental approximations from it. The new rule is 
semi-local in the cr-space, and, as a consequence, can be applied to all membrane topologies 
in a unified way, in marked contrast with previously known rules. Using our rule, one can 
construct functions corresponding to given matrices such that the fundamental approxima- 
tions hold well, provided that these functions exist. Whether these functions exist for given 
matrices can be also determined. As a physical application, for given matrices x a ,p a of the 
matrix model, one can construct the geometrical shape and the momentum densities of the 
membranes. 

The new rule includes the linear approximation (4.3), which is the analog of the Bohr- 
Sommerfeld condition. The linear approximation has lead us to the particular structure of 
the eigenvalue distribution, namely the branching phenomenon of the eigenvalue sequences. 
The eigenvalue sequences, which we have introduced in this paper, are subsets of the all 
eigenvalues whose members can be linearly approximated locally. From the analog of the 
Bohr-Sommerfeld condition, we have shown that the branching phenomenon reflects the 
information of topology. 

Thus, we have clarified the manner the information of topology manifests itself in the 
eigenvalue distribution. It is natural to further ask the question: "How completely can 
we read off the information of topology from given matrices?". We shall give here some 
observations which are essential to this question. In the first place, our argument implies that 
there is no information of topology in such ill-behaved matrices for which the fundamental 
approximation formulae (2.7)-(2.9) do not hold well. Indeed, it is only for the case (2.7)- 
(2.9) work, that the linear approximation (4.3) should hold. Hence, even the existence of the 
eigenvalue sequences is not guaranteed for those ill-behaved matrices. Secondly, there occurs 
overlapping of topologies when we consider the interaction of membranes. For a typical 
example, let us consider process shown in Fig. 8. At first there are two spheres. Then, these 
spheres approach each other and the distance A (in the target space) between two spheres 
reduces to zero gradually, and finally the two spheres merge into a sphere. This overlapping 
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Figure 8: Typical case of topology chang- 
ing of membrane 



Figure 9: Typical case of topology chang- 
ing from the viewpoint of eigenvalue dis- 
tribution 



is also present in the eigenvalue distribution. The process from the viewpoint of eigenvalue 
distribution is as follows. In Fig. 9 the eigenvalue distribution of the matrix corresponding 
to the height in Fig. 8 is shown. The eigenvalue distribution consists of two eigenvalue 
sequences. If the distance A between the two eigenvalue sequences gradually reduces to zero, 
then we cannot distinguish the eigenvalue distribution from that of a matrix corresponding 
to one sphere. It is interesting to treat topology changing processes of membranes by the 
matrix model in this way. 

Our discussion in this paper has been of purely kinematical nature. To explore the 
dynamical implication of our rule is also clearly important. For example, our consideration 
has made clear the distinction between the configurations of matrices which approximate 
membranes well and which do not. It is interesting to consider whether and how the former 
configurations dominate in the path integral of the matrix model. 

We would like to comment on the issue of the membrane instability [14] . Let us consider 
a configuration of membranes which has a spike-like portion whose area is less than 1/N. If 
we simply apply the analog of the Bohr-Sommerfeld condition, we should fail to include the 
information of the spike into the matrices. Stated more appropriately, our argument tells us 
that the configuration cannot be well approximated by N x N matrices. We want to stress 
that this spike has an essential difference to the spike which is considered in the membrane 
instability. One uses the word spike for a portion of a surface when its linear dimension is 
large, and at the same time its area is small. The difference between the spike in our context 
and the spike in the instability context lies in the meaning of the area. In the former, the 
area means area in the cr-space, that is essentially p + . In the latter, the area means area 
in the target space, or the energy of the spike. This difference is meaningful. Indeed, there 
are membranes which have portions that have small energy but large p + or vice versa. In 
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particular, one can construct configurations of the matrix model which approximate well the 
membranes with spikes in the sense of the membrane instability. 

As a direction of further investigation of the matrix regularization procedure itself, we 
recall the discussion at the end of section 4. Our rule in the present form does not include 
tunneling effects between sequences. Although our rule gives correct overall behaviour of 
the matrix elements, we could not trust the rule in the present form to investigate the 
matrix elements in the immediate vicinity of the merging and branching processes. Concrete 
examples of the processes are the points B, C in example 3 in section 5. 16 It is an important 
task to extend our rule to incorporate the tunneling effects[15]. One possible strategy would 
be to revisit the analog problem, namely the quantum mechanics of a particle in a double- 
well potential. We can construct a formula to relate the semi-classical wave functions in both 
wells, which is valid even for the energy level near the local maximum, extending the ordinary 
argument in the WKB approximation using Airy functions. Another interesting question is 
the uniqueness of the correspondence which gives the fundamental approximations. Although 
we cannot, at present, provide the proof, we suspect that the correspondence rule from which 
(2.7)-(2.9) can be derived is unique up to similarity transformation. Indeed, examples studied 
in section 5 suggest that the previously known rules and our new rule are the same up to 
similarity transformation. An immediate consequence of the uniqueness is that a change of 
the reference function should amount to a similarity transformation. 17 

It is believed that the matrix regularization can be extended for general even dimensional 
base spaces on which the Lie brackets can be defined. The correspondence rule between 
matrices and functions can be easily constructed by using tensor product of matrices, when 
the topology of the base space are given by direct product of some two dimensional spaces. 
These extensions are important in the matrix model of M-theory, in order to incorporate 
longitudinal 5-branes. To extend our analysis to study topological properties of these higher 
dimensional objects is also an interesting problem. We believe that the analysis analogous 
to the WKB approximation, used throughout in this paper, will also be useful for the higher 
dimensional case. However, it would be a challenging task, since the WKB approximation 
itself is not fully understood for generic non-integrable Hamiltonian systems on four or more 
dimensional phase spaces, compared to that for the necessarily integrable systems on two 
dimensional phase spaces considered in this paper. 

We conclude with three possible applications of the new correspondence rule. 

(1) One can construct various interesting configurations of the matrix model by our rule. 
A particular merit of our rule in this respect is that it can be applied to a membrane which 

16 We can trust our rule near appearing or disappearing processes, such as the point C in example 2 of 
section 5. The jump in the plot of the difference of eigenvalues is a natural consequence of our rule. 

17 This property implies that the fundamental approximations are no good for a configuration of matrices 
which consists of matrices corresponding to different topologies. 
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has genus higher than two, with no more difficulty than to a membrane having topology of 
a sphere or a torus. 

(2) Our rule may be useful when investigating the Lorentz symmetry of the matrix model. 
It has been tried to regularize the Lorentz generators of the continuum theory in order to 
construct those of the matrix model [16]. However, since some of the Lorentz generators are 
not built up by simple multiplications or integrations or Lie brackets, one has been inclined 
to use the basis expansion of the previously known correspondence rules. Since there are 
many different expansions for different topologies, it has been difficult to define the Lorentz 
generators in a unique way. Since our rule can be applied uniformly to all topologies, it is a 
promising tool in constructing definitions of the Lorentz generators of the matrix model in 
a unique way. 

(3) The geometrical interpretation of matrix elements in our rule may make the problem 
of the large N limit of the matrix model accessible. The problem can be interpreted as a 
renormalization of the membrane theory. We hope that our rule, by determining the short- 
distance (or rather, the small-area) degrees of freedom, enables us to construct a block-spin 
transformation of the matrix model. 

I would like to thank, first of all, Prof. T. Yoneya for discussions, encouragements and 
careful reading of the manuscript. I would like to thank, for discussions and encouragements, 
other members and former members of Komaba particle theory group, especially Y. Aisaka, 
S. Dobashi and T. Sato. Also, I would like to thank W. Taylor, S. Iso and T. Shimada for 
discussions and encouragements, and J. Hoppe for helpful and encouraging correspondence. 

Note Added: This paper is an extended version of the author's master thesis [17] sub- 
mitted to University of Tokyo on April 2002, where the basic results were preliminarily 
reported. The author has recently noticed that Hyakutake has constructed matrices corre- 
sponding to axial symmetric membrane configurations [18], which are special cases of the 
general prescription given in section 4 of this paper. 
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